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As computing resources are limited, choosing the parameters for a full Lattice QCD simulation 
always amounts to a compromise between the competing objectives of a lattice spacing as small, 
quarks as light, and a volume as large as possible. Aiming to push unquenched simulations with the 
Wilson action towards the computationally expensive regime of small quark masses we address the 
question whether one can possibly save computing time by extrapolating results from small lattices to 
the infinite volume, prior to the usual chiral and continuum extrapolations. In the present work the 
systematic volume dependence of simulated pion and nucleon masses is investigated and compared 
with a long-standing analytic formula by Liischer and with results from Chiral Perturbation Theory 
(ChPT). We analyze data from Hybrid Monte Carlo simulations with the standard (unimproved) 
two-flavor Wilson action at two different lattice spacings of a ~ 0.08 fm and 0.13 fm. The quark 
masses considered correspond to approximately 85 and 50% (at the smaller a) and 36% (at the 
larger a) of the strange quark mass. At each quark mass we study at least three different lattices 
with L/a =10 to 24 sites in the spatial directions {L =0.85-2.08 fm). 

PACS numbers: 11.15.Ha, 12.38.Gc 


I. INTRODUCTION 

It is in the nature of any numerical Lattice QCD calcu¬ 
lation that it can only be done at non-zero lattice spacing 
and in finite volume. Moreover, due to limited computing 
resources the typical quark masses currently employed 
are still substantially larger than the masses of the phys¬ 
ical quarks. In order to obtain physically meaningful 
predictions, extrapolations of lattice results to the con¬ 
tinuum, the infinite volume and to small quark masses 
are necessary. In the context of spectrum calculations 
one usually extrapolates in the lattice spacing and the 
quark mass, while the volume is preferably chosen such 
that its systematic effect on the masses can be largely ne¬ 
glected. The underlying assumption is that if the linear 
spatial extent L of a lattice with periodic boundary con¬ 
ditions is much larger than the Compton wavelength of 
the pion (e.g. if TOttL ^ 5, according to a common rule of 
thumb), then a single hadron H is practically unaffected 
by the finite volume (except that its momentum must be 
an integer multiple of 2'k/L). Its mass mH{L) in partic¬ 
ular will be close to the infinite-volume value defined at 
fixed lattice spacing and quark mass as 

rriH = lim mH{L). (1) 

L—^oo 

If the box size is decreased until the hadron barely fits 
into the box, the virtual pion cloud that surrounds the 
particle due to vacuum polarization is distorted, and pi- 
ons may be exchanged “around the world”. As a conse¬ 
quence the mass of the hadron receives corrections of or¬ 
der to its asymptotic value, which are small com¬ 

pared to the typical statistical errors as long as the lattice 
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remains sufficiently large. When L gets very close to the 
size of the region in which the valence quarks are con¬ 
fined, however, the the quark wave functions of the en¬ 
closed hadron are distorted and one observes rapidly in¬ 
creasing finite volume effects approximately proportional 
to some negative power of L. 

In the present work we explore the practical implica¬ 
tions of this picture by investigating, for various fixed val¬ 
ues of the gauge coupling and the quark mass, the actual 
volume dependence of simulated light hadron masses. 
Against the background of our ORAL project—whose 
name is an acronym for “Going Realistic And Light”— 
we ask in particular under which circumstances extrapo¬ 
lations in the lattice volume could be appropriate to ob¬ 
tain infinite-volume results from sub-asymptotic lattices, 
which would allow one to save valuable computing time. 
To this end we compare our data to various finite-size 
mass shift formulae available from the literature. 

While in past years the chiral extrapolation and the 
reduction of discretization errors have been at the center 
of many theoretical and numerical studies, there have, 
until recently, been rather few systematic investigations 
into the lattice size dependence of light hadron masses. 
These include, first of all, an analytic work of 1986 by 
Liischer [1] in which a universal formula for the asymp¬ 
totic volume dependence of stable particle masses in ar¬ 
bitrary massive quantum field theories is proven. Some 
years later Fukugita et al. [2-6] carried out a systematic 
investigation of finite-size effects in pion, rho and nucleon 
masses from quenched and unquenched simulations (with 
the staggered action). Related numerical studies with 
staggered quarks came also from the MILC collabora¬ 
tion [7-9]. Recently the systematic dependence of light 
hadron masses and decay constants on the lattice volume 
has been receiving renewed attention, see e.g. Refs. [10- 
22]. These studies include, on the one hand, a deter- 
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mination of the pion mass shift in finite volume using 
Liischer’s asymptotic formula with input from infinite- 
volume ChPT up to NNLO [11]. On the other hand, the 
finite-size mass shift of the nucleon has been calculated 
using relativistic Baryon ChPT in finite volume up to 
NNLO [10]. In the following Section II we will briefly 
summarize those results to which we will compare our 
numerical data in Section V. The details of the under¬ 
lying simulations and the determination of light hadron 
masses and other observables will be described in Sec¬ 
tions III and IV. 


II. FINITE-SIZE MASS SHIFT FORMULAE 

We consider a stable hadron H (= tt, N) on a four¬ 
dimensional hypercubic space-time lattice of spatial vol¬ 
ume and sufficiently large time extent T, with lat¬ 
tice spacing a set equal to unity for convenience. Both 
the bare coupling g and the quark mass held fixed, for 
large L the mass mniL) of the hadron is supposed to 
become a universal function of the product rriT^L in the 
finite-volume continuum limit (which is obtained by tak¬ 
ing g 0 and simultaneously L —>■ oo, while keeping 
niTrL constant). Since finite-size effects probe the system 
at large distances L they are insensitive to short-distance 
effects, so that this function should be independent of the 
form and magnitude of any ultraviolet cut-off. It is there¬ 
fore expected to hold also for finite lattice spacings. 

Attributing finite-volume effects at large L to vacuum 
polarization effects, Liischer’s formula [1] applied to QCD 
relates the asymptotic mass shift 

AmniL) = mH{L) - mn ( 2 ) 

to the (infinite-volume) elastic forward scattering ampli¬ 
tude Ft^h{v), where v is the crossing variable. For the 
pion it is given in terms of by [23] 

o pOO _ 

Am.(L) = dye-^^^^F^.{iy) 

+ 0(e-™^). (3) 

Because of m > a/3/2 the error term is exponentially 
suppressed compared to the first term. Due to the nega¬ 
tive intrinsic parity of the pion and parity conservation in 
QCD there is no 3-pion vertex, so that the term referring 
to a 3-particle coupling in the general formula of Ref. [1] 
is absent. At leading order in the chiral expansion the 
scattering amplitude is given by the constant expression 
Inserting this into Eq. (3) yields 

mTr{L) — ViiTr 

4(2^)3/2 /2 (m,L)3/2’ 1 1 

where Ki is a modified Bessel function, and the second 
expression follows from its asymptotic behavior, Ki{x) ~ 


3 Ki{mT^L) 

Stt^ /2 


j\fx^ for large x. In addition one can take existing 
NLO and NNLO chiral corrections to the infinite-volume 
amplitude F^^^^ into account and solve Eq. (3) numeri¬ 
cally [11]. We will consider the practical effects of such 
corrections more cfosely in Section V. 

Ref. [23] also quotes a finite-size mass shift formula for 
the nucleon that can be evaluated if the ttN scattering 
amplitude as known from experiment is inserted. In a 
sense this formula has been superseded, however, by a 
recent result derived from Baryon ChPT in finite volume 
by the QCDSF-UKQCD collaboration [10]. Using the 
infrared regularization scheme [24] they obtain 


Aa{L) 


ig\moml 

167r2/2 

X f dx ^ Aro^L|n|^mQa;2 -|- (1 — a;) ^ 

Jo n 

( 6 ) 


for the nucleon finite-size mass shift AmN{L) at 0{p^) in 
the p-expansion of the chiral Lagrangian. The constants 
gA and /tt are to be taken in the chiral limit, mo is the 
nucleon mass in the chiral limit and the pion mass m.^. 
parameterizes the quark mass via the Gell-Mann-Oakes- 
Renner relation. The pion decay constant is normal¬ 
ized such that its physical value is 92.4 MeV. Kq is a mod¬ 
ified Bessef function, and the sum extends over all spatial 
3-vectors n with integer components Ui, i = 1,2,3, ex¬ 
cept n = 0. Hi can be interpreted as the number of times 
the pion moves around the lattice in the i-th direction. 
At O(p^) an additional contribution to the mass shift 
AmAr(L) is given by 


A/L) = 


3m^ 

47r2/) 


■E' 


( 2 ci - C 3 ) 


-C2- 


Ki{\n\m^L) 

|n|m,rA 

K2{\n\m^Ly 


(|n|m^L)2 


( 7 ) 


where Ci, C 2 and C 3 are effective couplings and Ki and 
K 2 are again modified Bessel functions. The complete 
QCDSF-UKQCD result for the nucleon finite-size mass 
shift at NNLO reads 


mN{L) - niN = Aa{L) + Ab{L) + 0{p^). (8) 


To apply this formula to simulated lattice data, in 
Ref. [10] the parameters of the chiral expansion in ( 6 ) and 
(7) are taken partly from phenomenology and partly from 
a fit of numerical data for ttin from relatively fine and 
large lattices to the (infinite-volume) O(p^) formula [25] 


niN = mo — 4cim^ — 


e/(A)- 


3gi 

327r// 


m” 


9a 


327r2/2 

3g^ 

2567r/2TOg 


647r2/2 \^TOo 

— - 8 ci -I- C 2 -I- 4 c 3 
Too 


0 (r 




( 9 ) 
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where the counterterm e\ (A) is taken at the renormaliza¬ 
tion scale A. With all parameters fixed in this way, the 
formulae (6) and (7) provide parameter-free predictions 
of the finite-volume effects in the nucleon mass. Eq. (8) 
has already been demonstrated to work remarkably well 
for various volumes at pion masses of around 550 and 
700 MeV [10], and we will show in Section V that it is 
also capable of describing the volume dependence of our 
nucleon masses at pion masses from about 640 down to 
approximately 420 MeV. The QCDSF-UKQCD collabo¬ 
ration has shown that if the leading (jnj = 1) terms of 
their 0{p^) formula (8) are expressed in a form that cor¬ 
responds to Liischer’s approach [23], his nucleon formula 
is essentially recovered. A remaining numerical discrep¬ 
ancy has recently been identified as being due to a miss¬ 
ing factor of two in the so-called pole term of Liischer’s 
nucleon formula [17]. An important advantage of the 
formula (8) is that it is valid not just asymptotically but 
also at smaller values of L, because its subleading terms 
(Jnj > 1) account also for those virtual pions that cross 
the boundary of the lattice more than just once. 

Besides the formulae (3) and (8) we will confront our 
data also with the observation by Fukugita et al. that 
their pion, rho and nucleon masses from simulations with 
dynamical staggered quarks followed a power law, 

Am//(L) oc with n ~ 2..3, (10) 

rather than Liischer’s formula [4, 5]. Their result has 
been interpreted such that at smaller, sub-asymptotic 
volumes the leading finite-size effect originates from a 
distortion of the hadronic wave-function itself, contrary 
to the large-volume picture of a squeezed cloud of virtual 
pions surrounding a point-like hadron. 


III. SIMULATION 


A numerical investigation of finite-size effects in LQCD 
requires gluon field ensembles from several different lat¬ 
tice volumes at fixed gauge coupling and quark mass. In 
order to take benefit from our previous SESAM and TyL 
projects [26, 27] we have conducted supplementary simu¬ 
lations using again the standard Wilson action, the gauge 
part of which is given by the plaquette action 


X fiKu 


l--ReTrWl^\x) 


( 11 ) 


and the quark part by 

Sq = '^q{x)M{x,y)q{y), ( 12 ) 

x,y 


where 


M{x,y) 


— K 


i:['' - -ff,)Uf,{x)6^ + p,^y 


+ {1 + lf,)Ul{x - 


(13) 


We worked at two different values of the gauge cou¬ 
pling parameter, /? = 5.32144 (with k = 0.1665) and 
[3 = 5.6 (with K = 0.1575,0.158). The larger (3 cor¬ 
responds to the value used previously by SESAM/TyL. 
The smaller (3 and the corresponding k result from linear 
extrapolations of lines of constant mps/rny and mpgL in 
the (/3, K)-plane, based on SESAM/TyL data and aim¬ 
ing at mps/my < 0.5 and mpsL « 5 on a 16^-lattice 
[28]. For each of these /3,«;-combinations we have pro¬ 
duced unquenched gauge field configurations for at least 
three different lattice volumes with L varying between 
10 and 16, thus complementing ensembles from SESAM 
and TyL with L = 16 and 24, respectively. Generating 
the configurations periodic boundary conditions were im¬ 
posed in all four space-time directions for the gauge field, 
while for the pseudofermions we used periodic bound¬ 
ary conditions in the spatial directions and antiperiodic 
boundary conditions in the temporal direction. Beside 
the original SESAM TAG code that was used on a 512- 
processor APEIOO/Quadrics (QH4) we worked with an 
adapted version of the code on APEmille. There, a 128- 
processor crate was used to generate the 16^ x 32-lattices, 
while the 12^ x 32-lattices were produced on a unit of 32 
processors. On ALiCE [29], the 128-node “Alpha Linux 
Cluster Engine” at the University of Wuppertal, an opti¬ 
mized C/MPI-version (with core routines written in As¬ 
sembler) [30, 31] of the SESAM code ran on partitions 
of 16 (12^ X 32 and 14^ x 32 lattices) and 8 processors 
(10^ X 32 lattice). All the codes are implementations of 
the 4>-version [32] of the Hybrid Monte Carlo algorithm 
for two degenerate quark flavors. 

Tables I and II give an overview of the production runs 
we conducted within the GRAL project. For reference 
and convenience we also list the corresponding figures 
for previous SESAM/TyL runs here and in the follow¬ 
ing tables. The SESAM/TyL simulations at {(3,k,L) = 
(5.6,0.1575,16), (5.6,0.1575,24) and (5.6,0.158,24) have 
been used for our analysis of finite-size effects. Ex¬ 
cept for some early SESAM runs (or parts thereof) fea¬ 
turing an even/odd decomposition of the quark matrix 
M, in all simulations the BiCGStab algorithm was used 
with locally-lexicographic SSOR preconditioning [33, 34] 
for the inversion of the full matrix. The linear system 
{M^M)X = (j) was solved in two steps: first M^Y = (p 
was solved for Y and then MX = Y was solved for X. 
(iViter) in Table I denotes the average numbers of iter¬ 
ations the solver needed until convergence. Note that 
the slanted numbers quoted for runs on ALiCE refer 
to the solve of M^Y = p only, whereas in the case of 
runs on APE they refer to the full two-step solution of 
{M^M)X = (j). For a comparison a relative factor of ap¬ 
proximately two must therefore be taken into account. 
Like in the SESAM/TyL simulations the stopping accu¬ 
racy R = \\MX — (/jj/jjAj| was chosen to be i? = 10“® in 
all GRAL runs. 

Our APE programs additionally feature an implemen¬ 
tation of the chronological start vector guess proposed in 
Ref. [35]. In our GRAL simulations the depth of the ex- 
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TABLE I: The parameters of our simulations within the GRAL project together with those of previous SESAM and TxL runs 
which we have compiled here for reference (see also Table II). The chronological start-vector guess (csg) has been implemented 
only in our APE code. Note that in the case of ALiCE runs the numbers for (Mter) (slanted) refer to the solve of M^'Y = (p 
only, whereas in the case of APE runs they refer to the full two-step solution of (M^M)X = tj)- Note also that in order to boost 
the initially low acceptance rate of only 41% in the simulation at (/3, k, L) = (5.32144,0.1665, 16) we decreased both the step 
size 5t and the number of molecular dynamics steps at some intermediate stage of the simulation. 


p 

K 

lPt 

Preend. 

Nesg 

Nmd 

5t 

T 

P 

J- acc 

(Niter) 

(□) 

5.32144 

0.1665 

12^32 

11-SSOR 

- 

125 ± 20 

0.004 


0.71 

147(6) 

0.53949(14) 



14^32 

11-SSOR 

- 

125 ± 20 

0.004 

0.5 

0.64 

130(6) 

0.53879(15) 



16^32 

11-SSOR 

7 

200 ± 40 

0.005 


0.41 








125 ± 20 

0.004 

0.5 

0.65 

315(9) 

0.538290(65) 

5.5 

0.1580 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.77 

45(1) 

0.555471(45) 


0.1590 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.71 

85(1) 

0.558164(38) 


0.1596 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.61 

138(2) 

0.559745(58) 


0.1600 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.40 

216(3) 

0.560776(47) 

5.6 

0.1560 

16^32 

even/odd 

6 

100 ± 20 

0.010 

1.0 

0.82 

86(1) 

0.569879(25) 


0.1565 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.77 

90(1) 

0.570721(22) 


0.1570 

16^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.67 

133(1) 

0.571592(27) 


0.1575 

10^32 

11-SSOR 

- 

100 ± 20 

0.010 

1.0 

0.87 

63(1) 

0.573114(27) 



12^32 

11-SSOR 

7 

100 ± 20 

0.010 

1.0 

0.76 

146(2) 

0.572771(30) 



14^32 

11-SSOR 

- 

100 ± 20 

0.010 

1.0 

0.62 

79(1) 

0.572598(22) 



16^32 

even/odd 

11 

100 ± 20 

0.010 

1.0 

0.78 

293(6) 





11-SSOR 

3 

71 ± 12 

0.007 

0.5 

0.73 

160(6) 

0.572550(27) 



24 AO 

11-SSOR 

6 

125 ± 20 

0.004 

0.5 

0.80 

109(1) 

0.572476(13) 


0.1580 

12^32 

11-SSOR 

7 

125 ± 20 

0.008 


0.85 

150(5) 

0.573793(32) 



14^32 

11-SSOR 

- 

100 ± 20 

0.005 

0.5 

0.88 

113(1) 

0.573677(25) 



16^32 

11-SSOR 

7 

125 ± 20 

0.006 

0.75 

0.66 

302(5) 

0.573461(25) 



24 AO 

11-SSOR 

6 

125 ± 20 

0.004 

0.5 

0.62 

256^) 

0.573375(16) 


trapolation, A^csg, was not optimized, however, but rather 
fixed to 7. 

Both for decreasing quark mass and increasing lattice 
volume (all other parameters constant, respectively) we 
observe a drop in the acceptance rate as anticipated. 

Table II lists the total number of generated trajec¬ 
tories, IVtrajj the first IVtherm of which we attribute to 
the thermalization phase and therefore discard, so that 
we are left with N^qui equilibrium configurations, respec¬ 
tively. (In the thermalization phase of each production 
run we approached the respective target quark mass adi- 
abatically from larger quark masses. These initial tra¬ 
jectories are in general not counted here. An exception 
to this rule is the run at (/3, k, L) = (5.32144,0.1665,16) 
where a rather long initial tuning phase incorporated sev¬ 
eral changes of the simulation parameters.) A^conf con¬ 
figurations out of these, separated by AA^traj (up to a 
uniform, random variation) trajectories, have been ana¬ 
lyzed further. 


Autocorrelation 


A suitable estimator of the true autocorrelation (or 
autocovariance) function for a finite time-series At, t = 


1,... ,Tmc, is given by 


1 Tmc—^ 

s - A>J %+. - (TJ. 

S—1 

(14) 

where the use of the “left” and “right” mean-value esti¬ 
mators 


{^)l - 


1 


Tuc—i 


Tmc — t 




1 


Tmc — t 


TMC—t 

^r+t 


(15) 

in general leads to a faster convergence of C^{t) to the 
true autocorrelation function for Tmc ^ oo [36]. From 
fits of the estimator of the normalized autocorrelation 
function. 


p^(t) = C^(t)/C^(Q), (16) 


to an exponential we extract estimates for the exponen¬ 
tial autocorrelation time r^p, defined as 


exp 


lim sup 

t^OO 


t 

-log|p^(t)|’ 


(17) 


for A = □, Alter- Due to the notorious difficulty of de¬ 
termining autocorrelation times from short time-series we 
refrained from an elaborate optimization of the fit ranges, 
and the values given in Table III should be considered as 
rough estimates of the exponential autocorrelation times 
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TABLE II: Overview of run lengths, thermalization times, nnmbers of equilibrium configurations, numbers of analyzed con¬ 
figurations and the distances between them, machines used for the generation of the configurations, and the corresponding 
projects. (See text for more details.) 


p 

K 

L^T 

At,aj 

-^therm 

^equi 

AAt.aj 

^conf 

Machine 

Project 

5.32144 

0.1665 

12“32 

10100 

1500 

8600 

50 ±6 

170 

ALiCE 

GRAL 



14^32 

5900 

800 

5100 

40 ±4 

129 

ALiCE 

GRAL 



16^32 

15300 

8600 

6700 

40 ±4 

169 

APElOO/mille 

GRAL 

5.5 

0.1580 

16^32 

4000 

1000 

3000 

25 ±3 

119 

APEIOO 

SESAM 


0.1590 

16^32 

6000 

1000 

5000 

25 ±3 

200 

APEIOO 

SESAM 


0.1596 

16^32 

5500 

500 

5000 

25 ±3 

199 

APEIOO 

SESAM 


0.1600 

16^32 

5500 

500 

5000 

25 ±3 

200 

APEIOO 

SESAM 

5.6 

0.1560 

16^32 

5700 

600 

5100 

25 

198 

APEIOO 

SESAM 


0.1565 

16^32 

5900 

700 

5200 

24 

208 

APEIOO 

SESAM 


0.1570 

16^32 

6000 

1000 

5000 

25 

201 

APEIOO 

SESAM 


0.1575 

10^32 

16000 

2600 

13400 

48 ±4 

278 

ALiCE 

GRAL 



12^32 

8000 

700 

7300 

30 ±4 

243 

APEmille 

GRAL 



14^32 

8400 

1400 

7000 

30 ±4 

231 

ALiCE 

GRAL 



16^32 

6500 

1400 

5100 

25 

206 

APEIOO 

SESAM 



24® 40 

5100 

500 

4600 

25 

185 

APEIOO 

TyL 


0.1580 

12®32 

3000 

500 

2500 

24 ± 2 

103 

APEmille 

GRAL 



14®32 

9100 

1300 

7800 

40 ±4 

195 

ALiCE 

GRAL 



16®32 

6500 

1100 

5400 

30 ±4 

181 

APEmille 

GRAL 



24® 40 

4500 

700 

3800 

24 

158 

APEIOO 

TyL 


TABLE III: Measured autocorrelation times for the average 
number of solver iterations, Aiter, and the plaquette, □, for 
those (parts of) runs in which the fermion matrix was 11-SSOR 
preconditioned (see Table I). 


P 

K 

L^T 

-^iter 

' exp 

— -^iter 

' int 

' exp 

'int 

5.32144 

0.1665 

12®32 

52(8) 

58(5) 

53(9) 

50(3) 



14®32 

45(7) 

39(3) 

46(7) 

26(2) 



16®32 

97(14) 

105(7) 

77(9) 

54(4) 




94(22) 

77(12) 

74(10) 

57(5) 

5.5 

0.1580 

16®32 

21(5) 

20(1) 

27(3) 

16(1) 


0.1590 

16®32 

25(7) 

25(1) 

31(3) 

13(1) 


0.1596 

16®32 

74(12) 

38(2) 

37(11) 

17(1) 


0.1600 

16®32 

56(6) 

46(3) 

64(7) 

34(2) 

5.6 

0.1560 

16®32 






0.1565 

16®32 

29(6) 

19(4) 

7(1) 

5(1) 


0.1570 

16®32 

35(6) 

25(5) 

9(3) 

6(1) 


0.1575 

10®32 

62(7) 

28(2) 

15(3) 

5(1) 



12®32 

24(3) 

20(1) 

15(3) 

6(1) 



14® 32 

88(27) 

34(3) 

26(8) 

8(1) 



16®32 

47(7) 

33(4) 

18(6) 

7(4) 



24® 40 

51(7) 

36(4) 

11(2) 

7(3) 


0.1580 

12®32 

25(5) 

24(1) 

9(3) 

4(1) 



14® 32 

14(2) 

12(1) 

10(1) 

4(1) 



16®32 

43(15) 

24(2) 

14(4) 

8(2) 



24® 40 

61(19) 

50(5) 

20(10) 

20(2) 


only. We have checked, however, that the differencing 
method described in Ref. [36] yields consistent results. 
The measured values for r^p are generally larger than 
the integrated autocorrelation times that we measure 
as usual with the help of Sokal’s “windowing” procedure 
and which are also shown in the table. We use the finite 


sum 

-| Tout 

with a variable cut-off Tcut to estimate Plotting 

the resulting values against Tcut does, ideally, reveal a 
plateau for Tcut —*■ Tmc- If a plateau does not emerge, 
we typically either find a maximum, or r^^(Tcut) is 
monotonously rising. If there is a maximum, we choose 
the corresponding value as best estimate of Oth¬ 

erwise we reverse Sokal’s proposal to choose Tcut larger 
than 4 to 6 times rd^.: we assume to lie in the inter¬ 
val defined by the intersections of the straight lines with 
slopes Tcut/4 and Tcut/6, respectively, with the curve 

''’iiit(^cut)- 

Comparing autocorrelation times for runs with dif¬ 
ferent lattice volumes we find only a weak increase of 
the autocorrelation times with increasing volume. More 
striking is the difference in the autocorrelation times be¬ 
tween the simulations at (3 = 5.6 and /? = 5.32144. At 
the stronger coupling the relatively large autocorrelation 
times reflect the long-ranged statistical fluctuations that 
we observe in the corresponding time series. These fluc¬ 
tuations are more severe on the smaller lattices where, 
moreover, zero modes of the Dirac matrix start playing 
a role. On the largest volume at this (3 the situation is 
somewhat better: While the autocorrelation times are 
comparable to those on the smaller volumes, we see no 
indication of exceptional configurations on this (16^) lat¬ 
tice. 
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IV. PHYSICAL OBSERVABLES 


A. Static quark potential 


We calculated the static quark potential in order to 
determine the Sommer parameter [37] that we use to set 
the physical scale. For the SESAM/TyL runs at /3 = 5.5 
and (3 = 5.6 the Sommer radii Rq = r^/a as listed in Ta¬ 
ble IV have been published previously in Refs. [38] and 
[39], respectively. Since the lattice-size dependence of i?o 
is assumed to be small and as we want to have a common 
length scale for the different simulated lattice volumes at 
a given gauge coupling and quark mass, we have adopted 
the i?o-value from the largest available lattice, respec¬ 
tively, also for the smaller ones. 

In order to determine the Sommer radius for the 16^ 
lattice at (/?, k) = (5.32144, 0.1665) we measured the Wil¬ 
son loops W(R, t) with temporal extents of up to r = 8 
and spatial separations of up to i? = -\/3 • 7 « 12 lat¬ 
tice units on the same configurations that we used for 
spectroscopy. We employed the modified Bresenham al¬ 
gorithm of Ref. [40] to include all possible lattice vectors 
R for a given separation R = jRj. Using the spatial APE 
smearing as described in Ref. [41], we applied 

link ^ a X link -|- staples (19) 


to the gauge links of each configuration before calculating 
the Wilson loops. We used a = 2.3 and performed N = 
26 iterations, followed by a projection back into the gauge 
group [42]. 

The asymptotic behavior of the static potential V(R) 
for sufficiently large times r is given by 

WiR, t) ~ CiR) ( 20 ) 


so that one can define the effective potential 


Ueff(R,r) = In 


W{R,t) 
W{R,t + 1)' 


( 21 ) 


Fig. 1 shows the r-dependence of the effective poten¬ 
tial 14ff(i?, r) for various values of R. At r = 3 the 
effective potential is already largely independent of t 
while the statistical errors are moderate, so that we de¬ 
termined V (R) from a single exponential fit in the range 
[Tmin,Tmax] = [3,4]. As Can be Seen from Fig. 2, the re¬ 
sulting values for V(R) show the expected behavior. We 
observe no indication of string breaking [43] and therefore 
fit the data in the range [i?min, Rmax] = [2.5,8] to 


ViR) = Vo-^ + <jR. (22) 

The upper boundary of the fit range was set to i?max = 8 
because up to this value the data correspond nicely to 
the expected linear behavior, with small statistical er¬ 
rors. With Rmax held fixed the lower boundary was de¬ 
termined by investigating the i?min-dependence of rg for 
various values i?min ^ 2. (Below this value one observes 
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FIG. 1: The effective potential, 14ff(R, v), at (/3, k) = 
(5.32144, 0.1665) on the 16^ lattice, for selected values of R 
in the range r = 1,..., 4. Larger values of r are dominated 
by statistical noise. 



FIG. 2: The static quark potential obtained from Wilson 
loops at (/3, k) = (5.32144, 0.1665) on the 16® lattice. 


a violation of rotational symmetry due to the finite lat¬ 
tice spacing.) From the fit we obtained the following 
parameters for the potential (in lattice units): 

Vo = 0.8378(87), e = 0.440(18), a = 0.08187(97). 

The Sommer scale Rg, which is defined through the 
force between two static quarks at some intermediate dis¬ 
tance. 


R2 


dV 

dR 


R—Rq 


1.65, 


(23) 


was obtained from these parameters according to 


Rg — 


1.65- e 


(24) 
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Our result for the simulation at (/?, k, L) = 
(5.32144,0.1665,16) is given in Table IV. The quoted 
uncertainty corresponds to the statistical error. 

We used rg = 0.5 fm to set the scale in our simulations. 
The resulting physical values for the momentum cut-off 
a~^, the lattice spacing a and the box size L are dis¬ 
played in Table IV. The smallest simulated /? of 5.32144 
is associated with a lattice spacing of a = 0.13 fm, corre¬ 
sponding to a momentum cut-off of 1.52 GeV. While we 
have to watch out for potentially large 0{a) discretiza¬ 
tion errors at this coupling, the physical volume is the 
biggest of all simulated volumes. With a linear extension 
of slightly more than 2 fm it is comparable in size with 
the TyL lattice at (/3, k) = (5.6,0.1575). 

B. Hadron masses and decay amplitudes 

In order to extract meson masses and amplitudes 
we followed the standard procedure of computing zero- 


momentum 2 -point functions (x = (x, r)) 

(0t(r)0(0)) = ^(0t(a;)(!l(0)) (25) 

X 

for the following pseudoscalar and vector operators: 

P{x) = q'{x)-f 5 q{x), (26a) 

Ai{x) = q'{x)j5j4qix), (26b) 

Vkix) =q'{x)-fkq{x). (26c) 

For the nucleon we have used the octet baryon operator 

N{x) = eabc{qa{x)Cj5q'b{x))q''{x), (27) 


where a, 5, c are color indices and C = 7472 is the charge 
conjugation matrix. We employed the gauge invariant 
Wuppertal smearing [44] at the source only (Is) or at 
both source and sink (ss). In the previous SESAM and 
TyL simulations = 50 smearing steps were used with 
a weight of a = 4.0. These parameters were originally 
optimized for the 16^ SESAM lattice and then adopted 
for the larger 24^ TyL lattice, too. In order to adapt 
these parameters to our smaller lattices we have inves¬ 
tigated the effect of smearing on the various volumes. 
We applied the Wuppertal smearing procedure to point 
sources ^(°^(x) of size with L = 24,16,14,12,10. We 
set ^^°^(x) = 0 except for the point at (L/ 2 , L/ 2 , L/ 2 ), 
which we defined as the origin of the respective lattice 
and where we set (/)(°^( 0 ) = 1. Applying the smearing 
prescription N times to with all t7^(x) = 1 we plot 
the amplitude of the resulting wave function along 
the (arbitrarily chosen) ( 0 , 0 , l)-direction relative to its 
maximum at the origin, i.e. |(()(^)( 0 , 0 , a; 3 )p/|(()(^)( 0 )p, 
versus X 3 /L, for various values of N and a. On inspec¬ 
tion of the resulting wave function shapes we selected 
the parameters N and a for our simulated volumes so 
as to make the respective wave function profile look ap¬ 
proximately like the SESAM one. The chosen smearing 



FIG. 3: The parameters N, a for the Wuppertal smearing 
scheme were chosen such as to yield approximately the same 
wave function shapes for both the smaller lattices and the 
SESAM lattice {L — 16). 

parameters are listed in Table V, while the corresponding 
wave function profiles are displayed in Eig. 3. 

The masses and amplitudes of mesons were obtained 
by correlated least-y^ fits of the (time-symmetrized) cor¬ 
relators O = P,A 4 ,Vk, to the parameteri¬ 

zation 

f^{C, M) = C (e-^^ + , (28) 

where M = am is the mass and (7^1(0 |(!l|p)p/2M if 
Ip) is the zero-momentum state of the particle associated 
with O. In the case of the nucleon we fitted the cor¬ 
relator (anti-symmetrized in r and T — t) to the single 
exponential 

f^{C,M)=Ce-^\ (29) 

The optimal lower limits of the fit intervals, were 

found as usual by examining the y^/dof-behavior and 
the stability of the masses with respect to Tmin, and by 
inspection of the effective masses. The upper limit, Tmaxj 
was in general kept fixed at r/2 for mesons and r/2 — 1 
for the nucleon. 

The masses (in lattice units) of the pseudoscalar and 
vector mesons and of the nucleon are listed in Table VI. 
As the masses from the Is and the ss correlators are 
consistent we only quote the values extracted from the 
latter. The quoted errors are statistical in nature and 
have been estimated with the jackknife method (after 
suitable blocking of the data). Table VI also shows 
wps(A)T, the linear box size in units of the pseudoscalar 
correlation length l/mps(T), where tops(T) is the pion 
mass in the given finite volume. It should be borne in 
mind that for sub-asymptotic volumes this value is in 
general significantly different from mpsL, where mps is 
the pseudoscalar mass in infinite volume. At {(3, k) = 
(5.32144, 0.1665) we attain our lightest quark mass, with 






TABLE IV: Sommer scale and resulting momentum cut-off, lattice spacing and lattice size for ro = 0.5 fm. 


p 

K 

L^T 

ro/a 

a-i [GeV] 

a [fm] 

L [fm] 

5.32144 

0.1665 

16^32 

3.845(37) 

1.517(15) 

0.1300(13) 

2.081(21) 

5.5 

0.1580 

16^32 

4.027(24) 

1.5893(95) 

0.12416(74) 

1.987(12) 


0.1590 

16^32 

4.386(26) 

1.731(11) 

0.11400(68) 

1.824(11) 


0.1596 

16^32 

4.675(34) 

1.845(14) 

0.10695(78) 

1.711(13) 


0.1600 

16^32 

4.889(30) 

1.929(12) 

0.10227(63) 

1.636(10) 

5.6 

0.1560 

16^32 

5.104(29) 

2.014(12) 

0.09796(56) 

1.5674(89) 


0.1565 

16^32 

5.283(52) 

2.085(21) 

0.09464(93) 

1.514(15) 


0.1570 

16^32 

5.475(72) 

2.161(29) 

0.0913(12) 

1.461(20) 


0.1575 

16^32 

5.959(77) 

2.352(31) 

0.0839(11) 

1.343(18) 



24^40 

5.892(27) 

2.325(11) 

0.08486(39) 

2.0367(94) 


0.1580 

24^40 

6.230(60) 

2.459(24) 

0.08026(78) 

1.926(19) 


TABLE V: Smearing parameters. 


L 

10 

12 

14 

16 

24 

a 

1.0 

2.0 

3.0 

4.0 

4.0 

N 

20 

30 

40 

50 

50 


iTips/mv being close to 0.5. Using ro = 0.5 fm to set the 
physical scale the hadron masses of Table VI translate 
into the values listed in Table VII. This table shows also 
the physical box sizes L [fm] which have been calculated 
using the lattice spacing a from the largest available lat¬ 
tice, respectively (see also Table IV). The dimensionless 
quantity 


Mr = {rompsf (30) 


is another measure of the quark mass, since for niq —>■ 0 
the pion mass behaves like oc ^/mq. At the physical 
strange quark mass it gives Mr ~ 3.1 [45]. At those 
parameter sets where we have simulated several lattice 
volumes the value of Mr ranges between Mr « 2.65 and 
Mr ~ 1.13, corresponding to about 85 and 36% of the 
value for the strange quark mass. 

The (unrenormalized) pseudoscalar decay constant, 
which is defined on the lattice by (for p = 0) 

Z^(0lA4lPS) = MpsUps, (31) 


has been obtained from 





C ls 
A \ 


MpsCI 


(32) 


where we have used the fact that amplitudes for local 
source and sink {ll) can be obtained from Is and ss am¬ 
plitudes according to 


Cii 


~C^' 


(33) 


The pseudoscalar mass Mps = Mp in Eq. (32) has been 
taken from a fit of (P'I'(t)F(O))*^, the amplitudes 


(C^®) from a fit of the local-smeared (smeared-smeared) 

correlator (A 4 (t)A 4 ( 0 ))*^^*^\ 

In order to determine the (unrenormalized) quark mass 
Mq = aniq as defined via the PCAC relation on the lat¬ 
tice. 


MpsZa{0\A4\PS) 

2 Zp(OlPlPS)’ 

we used the relation 

7 M - Mil f^_Mll^ 

^ 2 ycl' 2 C^^\l CX 


(34) 


(35) 


where the renormalization constant is defined as Zq = 
ZpjZA and the pseudoscalar mass is taken to be the 
average 


Mps = i (Mj? -k M|,® -k Mjf -k Ml®). (36) 

Our results (in lattice units) for the unrenormalized 
pseudoscalar decay constant Pps/^A are displayed in Ta¬ 
ble VI. The normalization of the pseudoscalar decay con¬ 
stant is such that the physical value is = 92.4 MeV. 
The same table also shows the results for the bare PCAC 
quark mass ZqMq. 


V. VOLUME DEPENDENCE OF PION AND 
NUCLEON MASSES 

The three parameter sets at which we have data from 
several lattice volumes, namely (/?, k) = (5.6,0.1575), 
(5.6,0.158) and (5.32144,0.1665), are characterized by 
the quark mass, which in turn can be expressed in terms 
of the pion mass via the GMOR relation. We quote the 
pion mass measured on the largest lattice, respectively, 
when we refer to a particular simulation point (/3, k). We 
have investigated the volume dependence of the pion, the 
rho and the nucleon at pion masses (before continuum 
extrapolation) of approximately 643 MeV, 490 MeV and 
419 MeV in the ranges 0.85-2.04 fm, 0.96-1.93 fm and 
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TABLE VI: Masses of the pseudoscalar and vector mesons and of the nucleon, the pseudoscalar decay constant and the PCAC 
quark mass in lattice units. 


p 

n 

L^T 

nipsL 

Mps 

My 

Mjy 


ZqMq 

5.32144 

0.1665 

It^32 

3.18(8) 

0.2648(67) 

0.508(26) 

0.788(26) 

0.062(10) 

0.0106(33) 



14^32 

3.6(1) 

0.2577(87) 

0.518(12) 

0.779(16) 

0.0757(56) 

0.0152(18) 



16^32 

4.42(7) 

0.2760(42) 

0.4999(78) 

0.727(11) 

0.0843(62) 

0.0155(12) 

5.5 

0.1580 

16^32 

8.85(6) 

0.5534(39) 

0.6506(46) 

1.026(18) 

0.1073(42) 

0.0821(35) 


0.1590 

16^32 

7.09(4) 

0.4429(26) 

0.5529(54) 

0.8718(78) 

0.0945(23) 

0.0544(12) 


0.1596 

16^32 

5.89(4) 

0.3682(27) 

0.4902(52) 

0.7640(75) 

0.0815(16) 

0.03724(81) 


0.1600 

16^32 

4.89(5) 

0.3058(34) 

0.4547(61) 

0.703(10) 

0.0750(23) 

0.0279(15) 

5.6 

0.1560 

16^32 

7.15(4) 

0.4469(23) 

0.5365(36) 

0.8533(62) 

0.0843(19) 

0.0620(11) 


0.1565 

16^32 

6.32(6) 

0.3948(38) 

0.4989(54) 

0.785(10) 

0.0805(18) 

0.0467(15) 


0.1570 

16^32 

5.52(5) 

0.3452(29) 

0.4527(52) 

0.7095(90) 

0.0726(16) 

0.0391(15) 


0.1575 

10^32 

4.92(6) 

0.4919(55) 

0.587(20) 

1.042(20) 

0.0284(30) 

0.0209(32) 



12^32 

4.3(1) 

0.3576(89) 

0.494(12) 

0.817(16) 

0.0429(34) 

0.0249(19) 



14^32 

4.27(6) 

0.3048(44) 

0.4413(66) 

0.719(16) 

0.0566(30) 

0.0261(22) 



16^32 

4.49(6) 

0.2806(35) 

0.4036(68) 

0.6254(89) 

0.0626(26) 

0.0275(16) 



24^40 

6.64(6) 

0.2765(26) 

0.3944(38) 

0.5920(75) 

0.0646(18) 

0.02680(68) 


0.1580 

12^32 

4.6(1) 

0.387(12) 

0.535(17) 

0.882(25) 

0.022(12) 

0.0113(88) 



14^32 

4.13(8) 

0.2949(60) 

0.4677(90) 

0.717(19) 

0.0233(32) 

0.0099(17) 



16^32 

3.72(8) 

0.2325(51) 

0.371(13) 

0.622(12) 

0.0469(26) 

0.0141(11) 



24^40 

4.78(8) 

0.1991(33) 

0.3519(86) 

0.500(12) 

0.0602(39) 

0.0157(11) 


TABLE VII: Masses of the pseudoscalar meson, the vector meson and the nucleon in physical units (using ro = 0.5 fm). 


P 

K 

L^T 

L [fm] 

(romps)^ 

mps/my 

mps [GeV] 

my [GeV] 

mjv [GeV] 

5.32144 

0.1665 

12^32 

1.56(2) 

1.037(57) 

0.521(23) 

0.402(11) 

0.771(40) 

1.195(42) 



14^32 

1.82(2) 

0.982(68) 

0.497(20) 

0.391(14) 

0.786(20) 

1.182(26) 



16^32 

2.08(2) 

1.126(43) 

0.552(11) 

0.4188(75) 

0.759(14) 

1.104(20) 

5.5 

0.1580 

16^32 

1.99(1) 

4.97(20) 

0.8506(31) 

0.8795(81) 

1.0340(95) 

1.631(30) 


0.1590 

16^32 

1.82(1) 

3.77(12) 

0.8010(53) 

0.7666(64) 

0.957(11) 

1.509(16) 


0.1596 

16^32 

1.71(1) 

2.96(10) 

0.7512(51) 

0.6793(70) 

0.904(12) 

1.410(17) 


0.1600 

16^32 

1.64(1) 

2.235(85) 

0.6725(93) 

0.5901(75) 

0.877(13) 

1.356(21) 

5.6 

0.1560 

16^32 

1.567(9) 

5.20(18) 

0.8330(16) 

0.9002(69) 

1.0807(94) 

1.719(16) 


0.1565 

16^32 

1.51(1) 

4.35(25) 

0.7912(72) 

0.823(11) 

1.040(15) 

1.637(27) 


0.1570 

16^32 

1.46(2) 

3.57(21) 

0.7627(58) 

0.746(12) 

0.978(17) 

1.533(28) 


0.1575 

10^32 

0.849(4) 

8.40(59) 

0.838(30) 

1.144(14) 

1.365(48) 

2.424(47) 



12^32 

1.018(5) 

4.44(48) 

0.724(11) 

0.832(21) 

1.149(27) 

1.901(38) 



14^32 

1.188(5) 

3.22(18) 

0.691(11) 

0.709(11) 

1.026(16) 

1.671(39) 



16^32 

1.358(6) 

2.73(12) 

0.6952(99) 

0.6524(86) 

0.938(16) 

1.454(22) 



24^40 

2.037(9) 

2.654(90) 

0.7010(62) 

0.6429(67) 

0.9171(98) 

1.377(19) 


0.1580 

12^32 

0.963(9) 

5.80(88) 

0.722(20) 

0.951(30) 

1.316(44) 

2.167(66) 



14^32 

1.12(1) 

3.37(28) 

0.630(13) 

0.725(16) 

1.150(25) 

1.763(50) 



16^32 

1.28(1) 

2.10(15) 

0.627(21) 

0.572(14) 

0.912(33) 

1.530(32) 



24^40 

1.93(2) 

1.539(74) 

0.566(17) 

0.4896(94) 

0.865(23) 

1.228(31) 


1.56-2.08 fm, respectively. Due to angular momentum 
conservation the decay p ^ tttt is suppressed on small 
lattices where the minimum non-zero momentum 2tt jL 
is large. We therefore incorporate the rho resonance in 
our phenomenological analysis of finite-size effects, be¬ 
cause on the lattices considered here it should be stable. 


tive differences of the masses measured at L and Lmax, 
Mh{L) — Mh{L 

max) 


Rh{L) = 


MniL^ax) 


(37) 


where H = PS,V,N and Lmax = 24 (Lmax = 16) for 
P = 5.6 (/? = 5.32144). For both quark masses at /3 = 5.6 
we find a large variation of the hadron masses over the 


Figs. 4-6 show, for our three different quark masses, 
the pion, rho and nucleon masses in physical units as 
functions of the box-size. In Table VIII we list the rela¬ 


considered range of lattice sizes. While the finite-size ef¬ 
fects in the pion, rho and nucleon masses are relatively 
small (of the order of a few percent) if one compares 
only the two largest lattices at k = 0.1575, they rapidly 
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L[fm] 


L[fm] 


FIG. 4: Box-size dependence of the pseudoscalar and vec¬ 
tor meson masses and of the nucleon mass at (/3, k) = 
(5.6,0.1575). The solid lines result from fits to an exponen¬ 
tial, Eq. (38), while the dashed lines represent fits to a power 
law, Eq. (39). The curves are dotted outside the fit interval. 



FIG. 5: Fits as in Fig. 4 for (/3, k) = (5.6, 0.158). 


TABLE VIII: Ratios of the pseudoscalar mass and the chiral 
symmetry breaking scale, and the relative finite-size effects 
according to Eq. (37). 


d K 

L^T 

Mps 

Rps 

Rv 

Rn 

47rZ7^Eps 

5.32144 0.1665 

12-^32 

0.34(6) 

-0.04(3) 

0.02(5) 

0.08(4) 


14^32 

0.27(2) 

-0.07(3) 

0.04(3) 

0.07(3) 


16^32 

0.26(2) 

0 

0 

0 

5.6 0.1575 

10^32 

1.4(1) 

0.78(3) 

0.49(5) 

0.76(4) 


12^32 

0.66(6) 

0.29(3) 

0.25(3) 

0.38(3) 


14^32 

0.43(2) 

0.10(2) 

0.12(2) 

0.21(3) 


16^32 

0.36(2) 

0.01(2) 

0.02(2) 

0.06(2) 


24^40 

0.341(10) 

0 

0 

0 

0.1580 

12^32 

1.4(8) 

0.94(7) 

0.52(6) 

0.76(7) 


14^32 

1.0(1) 

0.48(4) 

0.33(4) 

0.44(5) 


16^32 

0.39(2) 

0.17(3) 

0.05(5) 

0.25(4) 


24^40 

0.26(2) 

0 

0 

0 


FIG. 6: Fits as in Fig. 4 for (/I, k) = (5.32144,0.1665). The 
fact that the pion and the rho show no monotonous increase of 
the finite-size shift towards decreasing box-size is ascribed to 
the smallness of the effect and statistics. (The simulations of 
the smaller lattices in particular at /3 = 5.32144 were affected 
by sizeable fluctuations.) In the case of the nucleon the finite- 
size effect is somewhat more significant. 


grow on the smaller volumes (~ 50-80% at L = 10). 
The rate of the increase is hadron dependent: while at 
large L the pion has the smallest relative finite-size ef¬ 
fect, the relative shift in the pion mass grows strongest 
with decreasing L, until it exceeds the effect in the rho 
mass from L = 12 and that in the nucleon mass from 
L = 10 downwards. Considering the finite-size effects at 
K = 0.158 (corresponding to a lower quark mass) we no¬ 
tice that at a given value of L the finite-size effects are 
generally much larger at k = 0.158 than at 0.1575. Again 
we observe that the pion is subject to the strongest rel¬ 
ative effect in the regime of small volumes. Finally, at 
{13, k) = (5.32144,0.1665) (corresponding to the lightest 
of our quark masses), we find rather small finite-size ef¬ 
fects of only a few percent in the simulated L-range, for 
all considered hadrons. In view of the small m-psL values 
(see Table VII) this is quite remarkable: if the finite-size 
effects would only depend on mpsT we would expect the 
effects at /3 = 5.32144 to be of about the same order of 
magnitude as those at the smaller volumes at (3 = 5.6. 
On the other hand, due to the large lattice spacing the 
volumes at /3 = 5.32144 are, in terms of the physical size, 
comparable with the larger volumes at (3 = 5.6. This 
strongly suggests that there are, in fact, different mecha¬ 
nisms responsible for the observed mass shifts, and that 
the transition between them is in our case characterized 
by the absolute physical lattice size rather than the prod¬ 
uct mpsL. Quite independently of the pion mass the re¬ 
gion in L where finite-size shifts start to become large is 
located at around 1.5 fm in our simulations. 
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A. Fits of the volume dependence 

First we attempted to describe the volume dependence 
of our simulated masses phenomenologically by fitting 
them to two different parameterizations. One of these 
parameterizations is inspired by Liischer’s exponential 
leading-order mass shift formula (5) for the pion, while 
the other one is directly given by the power law observed 
by Fukugita et al., Eq. (10). Although neither of these 
approaches can a priori be expected to be valid over the 
entire range of considered lattice volumes, and although 
Liischer’s formula strictly speaking has no free fit param¬ 
eters, on practical grounds it is still interesting whether 
based on either of the two functional forms an empir¬ 
ical description can be found that connects small and 
medium-sized volumes to the asymptotic regime. 

The curves in Figs. 4 and 5 show fits of the data for 
the pion, rho and nucleon masses {H = PS,V,N) to the 
exponential function 

mniL) = ran + cL~^ , (38) 

and, for comparison, to the power law 

mniL) = rriH + cL~^. (39) 

In the case of the pion (where mn = rnps) the mass 
Tops in Eq. (38) was treated as a fit parameter; the re¬ 
sult was used as input for the fits of the rho and the 
nucleon data, so that all the fits displayed in Figs. 4 and 
5 had two free parameters. As can be seen from the 
plots both parameterizations describe the data reason¬ 
ably well within the fit interval, but regarding the asymp¬ 
totic behavior the exponential ansatz is clearly superior. 
At (/?, k) = (5.6,0.1575), for example, all infinite-volume 
masses mn resulting from the exponential fits are com¬ 
patible with the data from the largest, 24^ lattice, which 
are assumed to bear no significant finite-size effects. In 
contrast, fitting the data to the power law yields numbers 
for rriH that grossly underestimate the true asymptotic 
masses. Varying the right boundary of the fit interval 
we find that for small box-sizes (L < 1.5 fm) where the 
finite-size effects are of the order of several percent the 
power law provides an acceptable description of the data. 
At the two larger quark masses this corresponds to the 
regime of mpsL < 4.5..4.8, in accordance with the com¬ 
mon rule of thumb that only for mpsL > 5 finite-size 
mass shifts are exponentially suppressed. (In the light 
of our results at (/?, k) = (5.32144,0.1665) it appears as 
if this rule of thumb could be relaxed as long as as the 
physical lattice extent L remains sufficiently large.) We 
find that as soon as we include data from larger volumes 
(where the mass shifts are small) into the fits the expo¬ 
nential ansatz yields better values for both y^/dof and 
rriH- In order to test whether this ansatz is suitable for 
an extrapolation from the small lattices to the infinite 
volume we fitted the masses at (/?, k) = (5.6,0.1575) and 
(5.6, 0.158) only up to L = 16. Assuming that the finite- 
size effects on the 24^ lattice are not significant it can be 



FIG. 7: Infinite-volume extrapolation of the masses at 
(/3, k) = (5.6, 0.1575). The solid lines correspond to exponen¬ 
tial fits according to Eq. (38). The curves are dotted outside 
the fit intervals. 



FIG. 8: Extrapolations as in Fig. 7 for {(3,k) = (5.6,0.158). 


seen from Figs. 7 and 8 that the asymptotic pion masses 
are generally underestimated considerably, while in the 
case of the nucleon mass the extrapolation works rather 
well. In either case the the extrapolation tends to yield 
lower bounds to the infinite-volume masses, the system¬ 
atic uncertainty of which can be estimated by varying 
the boundaries of the fit intervals. 

It should be mentioned that alternative fit formulae 
(obtained e.g. by changing the exponent of L from —3/2 
to — I in Eq. (38), introducing an additional variable fac¬ 
tor in the exponential of Eq. (38), or treating the expo¬ 
nent of L as a free parameter in Eq. (39)) may also be 
used to describe the data. They do not, however, lead to 
significant improvements and/or require even more free 
parameters. 

The main lesson we have learned from this exercise is 
the following: If one has got hadron masses from more 
than two different lattice volumes (at fixed coupling and 
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quark mass) and wants to estimate the infinite-volume 
masses on the basis of a fit, one should use an expo¬ 
nential ansatz rather than the power law. Extrapolating 
the exponential fit will produce lower bounds to the true 
asymptotic masses, and these bounds are generally better 
than those that can be obtained with the power law. 


B. Applicability of Chiral Perturbation Theory 


In order to better understand why it is problematic 
to extrapolate from small volumes to the infinite vol¬ 
ume on the basis of the simple formulae (38) and (39) 
one needs to appreciate their respective origin and scope. 
The power law (39) is supposed to originate from a dis¬ 
tortion of the hadron wave function (or from a modifica¬ 
tion of the effect of virtual particles traveling around the 
lattice by a model-dependent form factor that accounts 
for the finite hadron extent [4]) at quite small volumes. 
Consequently, the L“^-behavior is not expected to per¬ 
sist towards large volumes, which is in fact borne out by 
our data. On the other hand, the formula (38) essen¬ 
tially corresponds to Liischer’s asymptotic formulae for 
the pion. Liischer’s general formula for the volume de¬ 
pendence of stable particles in a finite volume represents 
the leading term of a large L expansion, meaning that 
whenever the relative suppression factor 

subleading ^ ^ 
leading V / 

is not small, subleading effects may be of practical rele¬ 
vance. 


1. Pion 


In the case of the pion we also rely on effective field 
theory to provide us with an analytic expression for the 
elastic forward scattering amplitude EVir- At leading or¬ 
der in the chiral expansion this amplitude is given by the 
constant expression = —rn^/f^. Inserting this into 
the Liischer formula eventually leads to Eq. (5), which 
has the functional form of our exponential ansatz (38). 
We have seen that the data can be described quite well by 
the parameterization (38). It is therefore instructive to 
compare our results for the parameter c to the constant 


C = 


m 


3/2 

PS 


4(2^)3/2 (Z;^Vps)2 


(41) 


(c/. Eq. (5)), where at each (/3, k) we take mps and 
from the largest available lattice, respectively 
(see Table VI). For our simulation points {(}, k) = 
(5.6,0.1575), (5.6,0.158) and (5.32144,0.1665) we have 
C/GeV^/^ = 1.089(61), 0.745(98) and 0.79(12), respec¬ 
tively. (Using = 137 MeV and = 92.4 MeV the 
natural value is C = 0.283 GeV“^^^.) Comparing the 


first two of these values to the results for c in Tables XI 
and XII (exp,PS) we see that the relative factor between 
c and C is 0(10) assuming that Za = 0(1). The dis¬ 
crepancy is generally larger for smaller values of the left 
fit boundary, Li, but decreases for increasing Li. 

The large differences between the coefficients c from 
the fits to our pion data and C from the Liischer for¬ 
mula (with LO ChPT input) reflect the fact that not 
all of our data sets for the different volumes at (/3, k) = 
(5.6,0.1575) and (5.6,0.158) comply with the conditions 
under which the application of this formula is justified. 
Recall that these conditions are: (i) sufficiently large lat¬ 
tice volumes (because the Liischer formula corresponds to 
the leading term in a large L expansion), and (ii) small 
pion masses (because we take the pion scattering ampli¬ 
tude from chiral perturbation theory). 

Quite recently, the finite-size shift of the pion mass 
has been determined using Liischer’s formula with the 
TTTT forward scattering amplitude taken from two-flavor 
chiral perturbation theory up to NNLO in the chiral ex¬ 
pansion [11]. These results have then been compared 
to the leading order chiral expression for the pion mass 
in finite volume (including the large-L suppressed terms 
neglected by the Liischer formula) in order to estimate 
the effect of subleading terms in the large L expansion. 
Both aspects of this investigation rely on chiral pertur¬ 
bation theory as an expansion in the pion mass ttItt and 
the particle momenta p, both of which have to be small 
compared to the chiral symmetry breaking scale that is 
usually identified with 47r/^. The conditions of applica¬ 
bility thus read 


and 


m-TT 

47r/^ 


<C 1 


(42) 


P 

47r/^ 


< 1 . 


(43) 


In a periodic finite box of size L, where the particles’ 
momenta can only take discrete values pfe = 27rnfe/L with 
nfe G Z, the second condition directly translates into a 
bound on the box size. 


A » — ~ Ifm, 


(44) 


where we have used the physical value of /,r- (As shown 
in Ref. [11] the pion mass dependence of /,r predicted by 
ChPT at NNLO is rather mild.) While a priori it is not 
clear what the practical significance of the relations (42) 
and (44) is, we can identify on the basis of Tables VIII 
and VII those data sets that stand the greatest chance 
of meeting these conditions. From Table VIII one can 
see that the ratios Mps(L)/(47rZ)(^Fps(L)) for all simu¬ 
lated volumes at (/3, k) = (5.32144,0.1665) are relatively 
small and compatible with each other. The correspond¬ 
ing ratio at (/?,«:) = (5.6,0.158) is also relatively small 
for L = Lniax (and, moreover, comparable to the num¬ 
bers at (5.32144,0.1665)), but the value for the second 
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FIG. 9: Volume dependence of our pion masses in the regime 
^ 1-3 fm. The circles, squares and diamonds represent 

our data at (/3, k) = (5.6,0.1575) {m„ — 643 MeV), {P,k) = 
(5.6,0.158) (m^ = 490 MeV) and (/?, k) = (5.32144,0.1665) 
(m,r = 419 MeV), respectively. The curves correspond to 
Liischer’s formula with input from ChPT at LO (dashed), 
NLO (long-dashed) and NNLO (solid). The dotted curves 
show the full LO chiral expression. The dash-dotted curve is 
the full LO result shifted by the difference between the NNLO 
and the LO Liischer formula. 


largest lattice is already significantly larger. Considering 
only the largest lattice, respectively, Mps/(47rZ(^^Tps) 
is largest at (/3, k) = (5.6,0.1575), but here the value at 
L = 16 is still consistent with the one at Lmax = 24. 
In view of this and recalling the relative finite-size mass 
shifts R (Table VIII) we infer from Table VII that for 
rriTr = 643 MeV and TOtt = 490 MeV we can trust ChPT 
at most on the largest volumes with L « 2 fm (and pos¬ 
sibly the 16^ lattice with L « 1.4 fm at = 643 MeV), 
while the lattices with L < 1.4 fm are most probably too 
small. At {(3,k) = (5.32144,0.1665), on the other hand, 
where m.^. = 419 MeV, all lattices are larger than 1.5 fm 
due to the relatively large lattice spacing, and hence ap¬ 
pear large enough for ChPT to be applicable. 

In order to corroborate these findings we checked how 
our simulated pion masses tops(T), for L > 1.3 fm, relate 
to the results of Ref. [11]. There, the chiral expression 
for the amplitude TVtt has been written as an expansion 
in powers of 

= 16^2 p ^2p^ p ^3p^ p ^ ( 45 ) 

where the parameter ^ is defined as 




/ 


2 


(46) 


Inserting the expansion (45) up to F 4 or Fq into Liischer’s 
formula (3) for the pion and using the chiral expression 
for the isospin invariant amplitude A of Ref. [46], the 
leading term in the large-L expansion is obtained up to 


NLO and NNLO in the chiral expansion. (Correspond¬ 
ingly, inserting (45) into (3) only up to F 2 yields the 
LO expression (4).) In order to calculate the predicted 
finite size shift for the pion numerically for our three 
different pion masses we need to know the respective nn- 
merical value of the expansion parameter In order 
to avoid the difficulties associated with the renormaliza¬ 
tion of the pion decay constant one can use the analytic 
expression for the pion mass dependence of which is 
known to NNLO in ChPT. If we take the pion mass from 
the largest lattice as a first approximation to the asymp¬ 
totic pion mass TOtt, respectively, we obtain the curves 
displayed in Fig. 9. The dashed curves correspond to 
Liischer’s formula (4) with F^^t^ from ChPT at leading 
order. The long-dashed and solid curves show the NLO 
and NNLO predictions, respectively. For comparison, the 
dotted curves show the full leading-order chiral expres¬ 
sion {Nf = 2) for the pion mass in finite volume, given 

by 


m.n(L) = 


00 

m{n) 

n—1 


Ki{y/nm-„L) 

\/nm-,^L 


(47) 


where the multiplicity m(n) counts the number of integer 
vectors n satisfying nf -b + n| = n [11]. Since the 
modified Bessel function Ki {x) falls off exponentially for 
large x, the sum in (47) is rapidly converging. For n = 1 
Eq. (47) corresponds precisely to the LO Liischer formula 
(4). Finally, the dash-dotted curve in Fig. 9 represents 
the best currently available estimate of the full finite-size 
effect, obtained by adding to the Liischer formula with 
Ftttt at NNLO the difference between Eq. (47) and the 
Liischer formula with Ft^tt at LO. 

The main conclusion we draw from the plot is that 
for all our three pion masses and for our lattices with 
L > 1.3 fm the finite-size effects predicted by ChPT are 
considerably smaller than our statistical errors. On the 
largest lattices with L ~ 2 fm the maximal predicted 
finite-size correction (corresponding to the dash-dotted 
curve in the plot) is about 0.3% for the lightest pion and 
0.05% for the heaviest one. This is in accordance with our 
presumption that for all practical purposes the finite-size 
effects in the pion masses are negligible on our largest 
lattices. At L ~ 1.3 fm the finite-size shift ranges be¬ 
tween 1% for the heaviest and about 3% for the lightest 
pion, which is of the order of the statistical uncertain¬ 
ties. For A > 1.3 fm the differences between the full 
one-loop ChPT result and Liischer’s formula with Ft^-k at 
LO are comparably small, indicating that here the use 
of Liischer’s asymptotic formula is indeed justified; the 
maximal difference in the relative effects is about 50% 
at A ~ 1.3 fm for the smallest pion mass. By contrast, 
the difference between the relative effects predicted by 
Luscher’s formula with Ft^-k at NNLO and LO amounts, 
for the same lattice size, to a factor of 3.2 for the lightest 
and 4.5 for the heaviest pion. 

Incidentally, a formula analogous to (47) exists also for 
the pion decay constant. (Recently also an asymptotic 
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FIG. 10: Nucleon mass data from various collaborations as a 
function of oc mg, including our data. The curve corre¬ 
sponds to a fit of the data represented by the open symbols to 
Eq. (9). These data points are from simulations on relatively 
large and fine lattices. Note that the fit result is consistent 
with the physical pion and nucleon masses (indicated by the 
star). 


formula d la Liischer has been derived for [12].) The 
only difference is that the relative hnite size effect is neg¬ 
ative and (for Nf = 2) four times larger than that of the 
pion mass: 


ML) = U 


1 - 4 ^ 


/ N Ki{y/nmTrL) 

2_^ m(n)- -= --— 

n—1 ^ 


(48) 


We have already seen that the volume dependence of our 
pion masses can be accounted for by chiral perturbation 
theory on the largest lattices at most, and there is no 
reason to believe that this should be different for the de¬ 
cay constant. But we can at least check whether we can 
recover the relative factor of minus four. Without go¬ 
ing into the details we just state here that while we find 
the finite-size effect of the pion decay constant indeed to 
be negative, its magnitude is, on the smaller lattices at 
(/3, k) = (5.6,0.1575) and (5.6,0.158), about the same as 
that of the pion mass shift; on the second largest vol¬ 
ume at {(3, k) = (5.6, 0.1575) the relative shift in /,r(T) is 
about twice as big as the shift in TO.n.(L). Unfortunately 
at (/3, k) = (5.32144, 0.1665), corresponding to our small¬ 
est quark mass, we can make no definite statement. 


2. Nucleon 

Regarding the nucleon mass, replacing the simple ex¬ 
ponential (38) by an ansatz corresponding more closely to 
Liischer’s nucleon mass shift formula of Ref. [23] might 
be considered as the natural next step towards a bet¬ 
ter description of the volume dependence. (Although, as 
we have seen, the ansatz (38) describes the data already 


quite well.) But since Liischer’s nucleon formula can be 
seen as a special case of the formula (7), let us instead 
confront our data for the nucleon mass directly with the 
formulae ( 6 ) and (7). Following Ref. [10] we fix gA and 
to the physical values gA = 1.267, = 92.4 MeV, 

and set the couplings C 2 and C 3 to C 2 = 3.2GeV“^, 
C 3 = —3.4GeV“^. The remaining parameters mg, ci 
and ei{X) (where the renormalization scale A is chosen 
to be 1 GeV) are taken from a Ht of data from various 
unquenched simulations with 

a < 0.15 fm, niTrL > 5 and < 800 MeV (49) 

to Eq. (9). In Ref. [10], data from the QCDSF [10], 
UKQGD [47], GP-PACS [48] and JLQGD [49] collabora¬ 
tions have been used. These data are plotted in Fig. 10 
(open symbols), complemented by the results from our 
largest lattices, namely the TyL results at (/3, k, L) = 
(5.6,0.1575,24), (5.6,0.158,24) and the GRAL result at 
(5.32144,0.1665,16). We also include the SESAM re¬ 
sult at (5.6,0.1575,16) in the plot, and recent results 
from CP-PAGS for small quark masses but from quite 
coarse lattices [50] (solid symbols). Although the con¬ 
ditions (49) are to some extent arbitrary we stick to 
them for definiteness. Gonsequently we refrain from re¬ 
peating the fits of Ref. [10] with our or the new GP- 
PAGS data, because only the TyL point at (/3, k, A) = 
(5.6,0.1575,24) meets all of the requirements in (49). 
Instead we quote the result of Ht 1 in Ref. [10] where 
Too = 0.89(6) GeV, ci = -0.93(5) GeV^ and e\{X = 
1 GeV) = 2.8(4) GeV“^, consistent with phenomenol¬ 
ogy. The corresponding curve is represented by the 
solid line in Pig. 10. The fact that the TyL point at 
(/?, K, A) = (5.6,0.1575,24) lies close to the curve with¬ 
out having been included into the fit hints to a small 
0{a) effect at this point. 

Note that we use the standard Wilson plaquette and 
quark action with errors at 0{a), whereas the data from 
the other collaborations have all been generated with 
0 (a)-improved actions. 

The other TyL point at (/?, k, A) = (5.6,0.158, 24), cor¬ 
responding to a smaller pion mass, lies somewhat below 
the curve. Gorrecting it for the presumed finite-size ef¬ 
fects in the pion and the nucleon mass would shift it 
even slightly further away from the curve (recall that 
in this regime of larger A the finite-size effect is bigger 
for the nucleon than for the pion). The SESAM point 
at (/?, K, A) = (5.6,0.1575,16) illustrates how finite-size 
effects appear in such a plot. Correcting it for the finite- 
size effects in the pion and the nucleon masses (see Ta¬ 
ble VIII) would shift it to the lower left, towards the 
corresponding TyL point with A = 24. The data points 
from our largest lattices generally tend to lie somewhat 
below the curve, and this is also true for the GRAL point 
with {P,k) = (5.32144,0.1665). In view of the fluctua¬ 
tions in mTr{L) at this parameter set we plot in Fig. 10 
the mean of the respective pion masses at A = 12,14,16, 
with a corresponding error bar along the axis. Even 
with this uncertainty taken into account the deviation of 
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FIG. 11: Volume dependence of the nucleon mass for rtiT: = 
643 MeV. The dashed curve represents the 0(p^) term only, 
while the solid curve also includes the 0{p^) contribution. 
The dotted curve results if the pion mass is reduced by 10% 
in the O(p^) formula. 



FIG. 12: Same as Fig. 11 for rriT^ — 490 MeV. 


the GRAL point from the fit curve is significant. Consid¬ 
ering the relatively low cut-off of only about 1.5 GeV at 
this point (to be compared to a nucleon mass of 1.1 GeV) 
discretization errors might be responsible for the devia¬ 
tion. In case of the TyL data cut-off effects are expected 
to be less important, due to the smaller lattice spacings 
in these simulations. 

Using the parameters corresponding to the solid curve 
in Fig. 10 we can evaluate the finite-size formulae (6) 
and (7) and compare the results to our data. Our three 
sets of simulations with different lattice sizes correspond 
to pion masses of approximately 643, 490 and 419 MeV. 
Note that the latter two masses are lighter than the light¬ 
est of the pion masses investigated in Ref. [10] (732, 717 
and 545 MeV). The curves in Figs. 11-13 have been com¬ 
puted from equations (6) and (7) with no free parameters. 
Like in Ref. [10] the solid curves correspond to the 0(p‘*) 


FIG. 13: Same as Fig. 11 for rrin ~ 419 MeV. 
prediction 

mN{L) = niN +/^a{L) + Ab{L), (50) 

where ttim has been determined such that the calcu¬ 
lated value mjv(Lmax) equals the simulated mass from 
the largest lattice with L = respectively. Corre¬ 

spondingly, for the pion masses niTr we also take the sim¬ 
ulated value from the largest lattices. For the dashed 
curve, corresponding to the 0{p^) prediction, the O(p^) 
contribution from Af, in (50) has been omitted, while 
TOjv has been left unchanged. For all our pion masses we 
find a surprisingly good overall description of our data by 
the O(p^) prediction even down to lattice sizes of about 
1 fm. Replacing m,r from the largest, L = 16 lattice at 
(/?, k) = (5.32144,0.1665) by the mean of the pion masses 
from the L = 12,14,16 lattices (as we did in Fig. 10) does 
not lead to a significant difference in the resulting curve. 
Since both the statistical and the theoretical errors of the 
simulated TOAr(L) are smallest for the largest lattice, we 
consider the 0{p^) finite-size corrected nucleon mass 

7TiAr(L) = mAr(L) - Aa(L) - A{,(L), (51) 

taken at L = Lniax, as our best estimate of the asymptotic 
nucleon mass. Table IX shows the predicted infinite- 
volume masses for our simulations. The last column 
gives the relative mass shift on the largest lattice, re¬ 
spectively. Just as it was the case for the pion, the finite 
size effect in the nucleon at (/?, k, L) = (5.6,0.1575,24) 
is considerably smaller than the statistical uncertainty. 
At {P,k,L) = (5.6,0.158,24) and (5.32144,0.1665,16), 
on the other hand, the finite-size effects according to 
Eq. (51) amount to about 2% of the respective asymp¬ 
totic mass, which is comparable to the statistical errors. 

Gompared to a fit-based extrapolation the advantage 
of a formula without free parameters is of course that 
one can directly calculate the amount by which one has 
to shift the nucleon mass in order to compensate for the 
finite-size effect associated with a given volume, and that 







16 


TABLE IX: Finite-size shifts of the nucleon masses on the largest lattices as inferred from Eq. (50). A = {niN — 'ifiN)/mN is 
the relative deviation of the Monte Carlo value rriN from the shifted mass rfiN = ruN — Aa — A6 (all values to be taken at 
Lmax). We consider mjv(Lmax) as the best estimate of the true asymptotic mass. 


p 

n 

uips [GeV] 

^iv(-^/max) [GbV] 

^iV (-i'max) [06^] 

A 

5.6 

0.1575 

0.6429(67) 

1.370(19) 

1.377(19) 

0.53% 


0.1580 

0.4896(94) 

1.204(31) 

1.228(31) 


5.32144 

0.1665 

0.4188(75) 

1.081(20) 

1.104(20) 



one has control over the error. In practice, however, a 
remaining caveat is that the infinite-volume pion mass 
must be known. If one is working in a parameter regime 
where the finite-size effect in the pion mass is small (of 
the order of a few percent) one can apply the results of 
Ref. [11] to obtain an estimate of the true asymptotic 
mass. If this is unclear, but data from several (more 
than two) different volumes are available, one might still 
revert to an exponential fit and extrapolate. Since we 
have seen that such a “naive” extrapolation systemati¬ 
cally underestimates the true infinite-volume pion mass 
we illustrate, as an example, the impact of a by 10% 
smaller pion mass by the dotted curves in Figures 11-13. 
Although relative to the very nucleon mass shift the sys¬ 
tematic error associated with the uncertainty in the pion 
mass grows with L, its absolute value becomes less and 
less significant compared to the statistical errors of the 
data. On the one hand this means that (assuming the 
formula to exactly reproduce the volume dependence of 
the data and the statistical uncertainties to be all of com¬ 
parable size) in order to predict the asymptotic nucleon 
mass correctly (within the statistical errors) one needs 
to know the more accurately the smaller the physical 
size of the largest available lattice. If, on the other hand, 
L is sufficiently large so that one can reliably extrapolate 
the pion mass, the asymptotic nucleon mass can be de¬ 
termined quite accurately, already on the basis of a single 
lattice. 

C. Spatial Polyakov-type loops 

At our two larger quark masses we have observed in all 
considered quantities (pion, rho and nucleon masses, pion 
decay constant) a drastic increase of the finite-size shifts 
below a lattice size of approximately 1.5 fm. But above 
this size the finite-size effects are relatively small, also 
at our smallest quark mass. As we will now show, this 
kind of transition behavior is reflected in the behavior of 
spatial Polyakov-type loops. 

In the absence of dynamical quarks, i.e. in quenched 
QCD, the expectation value of the Polyakov loop, which 
is defined as 

(^)=(]^ET>'nC^4(x,r)^, (52) 

is zero in the confined phase, while in the deconfined 
phase (P) 7 ^ 0. Therefore in pure gauge theory (P) 


TABLE X: Ensemble averages of the real and imaginary parts 
of the mean spatial Polyakov loop (in 2 -direction). 


P K 

L^T 

(RePQ 

{Ira Pz) 

5.32144 0.1665 

12^32 

-0.00056433(3) 

-0.0003573(2) 


14®32 

-0.0000179(2) 

-0.00003209(4) 


16^32 

0.00012251(3) 

-0.00021856(4) 

5.6 0.1575 

10®32 

-0.0131032(3) 

-0.006341(8) 


12^32 

-0.0029148(3) 

0.0012405(6) 


14^32 

-0.00083854(3) 

0.00008748(3) 


16^32 

0.00020948(3) 

0.00004254(5) 


24^40 

-0.000078415(8) 

-0.00009462(1) 

0.1580 

12®32 

-0.003798(2) 

0.005014(3) 


14®32 

-0.00135725(6) 

0.0003214(3) 


16® 32 

-0.00037040(4) 

-0.00022400(4) 


24® 40 

-0.00010810(2) 

0.00012503(1) 


is an order parameter for the deconfining phase tran¬ 
sition. This is due to the global Z{3) symmetry of the 
pure SU (3) gauge theory which is spontaneously broken 
at the phase transition. In full QCD the Polyakov loop 
is not an order parameter because the Z{3) symmetry 
of the gluonic action is explicitly broken by the quark 
action, so that (P) is not exactly zero in the hadronic 
phase. In our simulations it is not the time extent T 
which is varied, but the spatial lattice size L. Due to the 
space-time symmetry of the Euclidean metric, however, 
similar considerations also apply to Polyakov-type loops 
in the spatial directions (also known as Wilson lines). 
Let us, for definiteness, consider the mean Wilson line in 
2 -direction, which is defined configuration-wise as 

1 ^ 

Pz = J^ E Yi U3ixi,X2,X3,T). (53) 

X1,X2,T = l 

As can be seen from Table X, the expectation values 
(Pz) for all GRAL simulations are indeed significantly 
different from zero, even on the largest lattices. While 
for the larger lattices the deviation from zero is relatively 
small it becomes more pronounced as the lattices shrink. 
(RePz) in particular takes increasingly negative values 
towards the smaller volumes. This can be understood by 
looking at the distribution of Pz- 

Fig. 14 shows the distribution in the complex plane 
of Pz for the lattice volumes simulated at (/?, k) = 
(5.6,0.1575). The lines in the plots indicate the three 
Z{3) directions 1, and Apart from the L- 

dependent fluctuations we observe for the larger lattices 
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FIG. 14: Distribution in the complex plane of the spa¬ 
tial Polyakov loop for the different lattices simulated at 
{/3,k) = (5.6,0.1575). The lines indicate the three Z{3) di¬ 
rections 1, and 



L[fml 


FIG. 15: Absolute values of the complex ensemble averages 
of the spatial Polyakov loop Pz. 


an approximately point-symmetric accumulation of the 
Wilson line around zero. This is reflected in the small¬ 
ness of I (Pz) I and the corresponding statistical errors at 
large L (Table X). The situation is somewhat different for 
the smallest, 10^ lattice, where the distribution of Wilson 
lines is visibly shifted towards the Z(3) directions 
and which leads to the relatively large negative 

value of (ReP 2 ). 

This shift can be understood e.g. by means of the 3-d 
Potts model with magnetic field that is recovered when 
the full QCD action is expanded first in the gauge cou¬ 
pling (3 and then in inverse powers of the sea quark mass 
[51]. Introducing the quark action in QCD is then equiva¬ 
lent to switching on a magnetic field h in the Potts model 
that breaks the Z (3) symmetry of the system. Consider¬ 
ing the phase of the spatial Polyakov loop as a spin that 
can take one of the three possible values 1, and 

g47ri/3^ the magnetic field aligns the spins to preferred 
directions depending on the sign of h: for h = — \h\ (cor¬ 
responding to antiperiodic spatial boundary conditions 
for sea quarks) the positive real axis is favored, whereas 
for h = -l-|ft.| (periodic spatial boundary conditions for 
sea quarks) the two directions and (point¬ 

ing towards negative values) are preferred. If we recall 
that we have used periodic spatial boundary conditions 
for sea quarks this explains the plot for L = 10 in Fig. 14. 

The implications of such a shift for finite-size effects 
in hadron masses have been explained in detail e.g. by 
Aoki et al. in the context of their comparative study of 
finite-size effects in quenched and full QCD simulations 
[2]. We briefly recapitulate their argument for our choice 
of boundary conditions: Let us consider a meson propa¬ 
gator r(r) on a lattice of size with a sufficiently large 
time extent T. A hopping parameter expansion of r(T) 
yields a representation of the meson propagator in terms 
of closed valence quark loops C going through the me¬ 
son source and sink. If we denote the corresponding link 


factors Trf][;g(^ Ui by P{C) for Polyakov-type loops that 
wind around the lattice in a spatial direction, and W (C) 
for ordinary Wilson-type loops, the meson propagator 
can be written as 

-(r(r)> = E «vT(VF(C)) + E Vvai(P(C)), 
c c 

(54) 

where L{C) is the length of the respective loop and the 
sign factor aval is equal to -1-1 for the periodic spatial 
boundary conditions used for valence quarks in our sim¬ 
ulations. From the discussion above we know that in 
the case of periodic spatial boundary conditions for both 
sea and valence quarks the contribution of Polyakov-type 
loops to the meson propagator (54) is negative. Since 
mean values of Wilson-type loops are always positive, 
the two contributions in (54) have opposite sign, which 
leads to a faster decrease of the correlator and thus to a 
larger meson mass. 

For fixed sea and valence quark mass this effect grows 
weaker for increasing lattice size because the contribution 
of the Polyakov-type loops decreases. This can clearly 
be seen from Fig. 15 for our larger pion masses. On 
the other hand, in a fixed lattice volume with periodic 
boundary conditions finite-size effects in hadron masses 
get increasingly significant both for decreasing sea and 
valence quark mass. This has been observed e.g. in 
Ref. [52], where partially quenched chiral extrapolations 
of the pseudoscalar and vector masses were studied for 
the various sea and valence quark masses at /3 = 5.5 and 
5.6. The sea quark mass dependence of the expectation 
value of the Wilson line can also be seen in Table X if 
one compares at /3 = 5.6 the value for a given lattice size 
at K = 0.1575 with the corresponding value at k = 0.158. 
We find that in the same volume (Re Pz) is more negative 
for the larger k, corresponding to a smaller quark mass. 
This leads to a stronger cancellation of the two terms in 
(54) and hence to a larger finite-size effect in the hadron 
masses. 
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FIG. 16: Same as Fig. 14 for (/3, k) = (5.6,0.158). 


-0.04 -0.02 0 0.02 0.04 -0.04 -0.02 0 0.02 0.04 -0.04 -0.02 0 0.02 0.04 

0.04 

0.02 

0 

-0.02 

-0.04 

FIG. 17: Same as Fig. 14 for (/I, k) = (5.32144,0.1665). 



For completeness. Figs. 16 and 17 show the distribution 
of the Wilson line for the simulated volumes at (/?, k) = 
(5.6,0.158) and (5.32144,0.1665). 

VI. DISCRETIZATION ERRORS 

One method that lends itself to checking on the signifi¬ 
cance of lattice artefacts in our simulations is to consider 
the PCAC relation 

i9^A“(a:) = 2mP°‘{x) (55) 

between the isovector axial current 

=q{x)-ff,j5^T°'q{x) (56) 

and the associated density 

P“(x) = g(a:)75ir“g(x), (57) 

where t“ denotes a Pauli matrix acting on the flavor in¬ 
dices of the quark field q. On the lattice, the bare quark 
mass ZqMq can be extracted from ratios of correlation 
functions, 


ZqMq 


1 {d^A-{x)J-) 

2 (P“(x)J“) 


0{a), 


(58) 


where the (smeared) source J“ is a suitable polynomial in 
the quark and gluon fields, and Zq = Zp/Z a- The PCAC 



FIG. 18: Box-size dependence of the relative shift in the 
PGAG quark mass at (/!,«;) = (5.6,0.1575). 
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FIG. 19: Same as Figure 18 for {13, k) = (5.6,0.158). 


relation (55) is an operator identity that holds for the 
Wilson action up to 0{a) effects. Consequently, its lat¬ 
tice version holds—up to those effects—for any choice of 
boundary conditions, source operators and lattice sizes. 
This means in particular that at fixed f3 and k any resid¬ 
ual lattice size dependence of the PCAC quark mass (58) 
must be a lattice artefact. Figs. 18-20 show the volume 
dependence of the relative deviation of the PCAC quark 
mass niq (L) from its value mq on the largest lattices, for 
our three (/3, k) combinations. 

We find that at {P,k) = (5.6,0.1575) the discretiza¬ 
tion errors appear to be small for L > 1 fm, while 
for {P,k) = (5.6,0.158) and (5.32144,0.1665) they are 
small only for L > 1.3 fm and L > 1.8 fm, respec¬ 
tively. On the smaller lattices the cut-off shows up in 
quark mass shifts of 20-40% (with large error bars on 
the smallest lattices). The fact that the cut-off effects 
are small for (/?,«:) = (5.6,0.1575) and somewhat larger 
for {(3,k) = (5.6,0.158) and (5.32144,0.1665) is consis¬ 
tent with our observations in Section V, where in Fig. 10 
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FIG. 20: Same as Figure 18 for (/9, k) = (5.32144,0.1665). 


we saw no significant lattice artefacts in the nucleon 
mass for m,r = 643 MeV, while for m,r = 490 MeV 
and niTr = 419 MeV the nucleon mass displayed some 
deviation from the curve (which represents a fit to 0{a)- 
improved data). 


VII. SUMMARY 

In order to investigate finite-size effects in stable light 
hadron masses obtained from Lattice QCD with two 
dynamical flavors of Wilson fermions we have com¬ 
plemented previous SESAM/TyL simulations at quark 
masses corresponding roughly to 85 and 50% of the 
strange quark mass by several runs at the same coupling 
of /3 = 5.6, but on smaller lattices. In addition we have 
carried out exploratory simulations using three different 
lattice volumes at a stronger coupling of /3 = 5.32144 
in an attempt to push our analysis towards the regime 
of lighter quark masses. We have succeeded in simulat¬ 
ing near nis/i, which in terms of rriT^/mp is close to the 
rho decay threshold of 0.5. The physical extent of the 
investigated lattices ranges between 0.85 fm and 2.04 fm. 

We have addressed, from a practical point of view, the 
question to what extent the volume dependence of the 
computed pion, rho and nucleon masses can be param¬ 
eterized by simple functions, and if with these functions 
an extrapolation from small and intermediate lattices to 
the infinite volume is possible. To this end we have 
compared an exponential ansatz motivated by Liischer’s 
mass-shift formula for the pion to the power law observed 
by Fukugita et al.. On the basis of various fits we con¬ 
clude that while the power law may be used to describe 
the volume dependence of the masses at volumes smaller 
than roughly (1.5 fm)^, over the full range of simulated 
lattices—and in particular with respect to the asymptotic 
behavior—the exponential ansatz is more appropriate. 
Although the extrapolation of simple exponential fits to 
the infinite volume in general provides only a lower bound 


to the asymptotic mass, this bound may be close to the 
true asymptotic value if the relative difference between 
the masses from the largest volumes incorporated in the 
fit is already quite small (of the order of a few percent). 
For small volumes alone, however, this is in general not 
the case. 

The exponential parameterization corresponds in its 
functional form precisely to Liischer’s asymptotic formula 
for the pion (with input from infinite-volume ChPT at 
leading order). Although we have found that the single 
exponential allows for a reasonable phenomenological de¬ 
scription of our light hadron masses over a wide range of 
lattice volumes, a large coefficient multiplying the expo¬ 
nential attests to the fact that the data points from the 
small lattices lie outside the parameter regime in which 
the original formula holds. In the case of the pion we have 
illustrated this by a comparison of our data to Liischer’s 
formula with input from ChPT up to NNLO and to the 
full LO ChPT result for the pion mass in finite volume. 
Of course, if an appropriate analytic prediction with con¬ 
trolled errors is available it is preferable to an extrapola¬ 
tion based on a fit with free parameters. We have found, 
however, that in the parameter regime of our simulations 
even the best currently available estimate for the pion 
finite-size effect, based on a combination of the asymp¬ 
totic Liischer formula with NNLO ChPT input and the 
full finite-volume (but LO) ChPT result, predicts mass 
shifts of a few percent only. This is comparable in size to 
the typical statistical errors and therefore hard to detect 
in practice. Our simulations at /? = 5.32144 probably 
are in a pion mass regime where the box-size dependence 
can be described by such a formula, but more statistics is 
needed to corroborate this assumption. While Liischer’s 
formula with input at next-to-leading and next-to-next- 
to-leading order ChPT can be used to control the conver¬ 
gence of the chiral expansion, a full higher order result 
from ChPT for the pion mass in finite volume would be 
highly useful to fully assess the role of the sub-leading 
terms in the large-L expansion. 

For the nucleon, a promising finite-size mass formula 
has recently become available from relativistic Baryon 
ChPT. We have shown for three different pion masses 
(two of which are smaller than the ones considered by 
the QCDSF-UKQCD collaboration) that it describes our 
simulated nucleon masses remarkably well even down to 
box-sizes of about 1 fm. We have also seen that above this 
size it can, in principle, be used to estimate the infinite- 
volume mass already on the basis of a single measure¬ 
ment, provided that the corresponding asymptotic pion 
mass is known. If, as in our case, data from several lattice 
volumes are available, they can be combined to obtain a 
reliable estimate with controllable errors. 
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TABLE XI 

Fit parameters for (/3, k) 

= (5.6,0.1575). 

d = 2 for “pow” 

(Eq. (39)) and d = 

1/2 for “exp” 

(Eq. (38)). 

Fit type 

H 

[Li,L2\ 

niH [GeV] 

c [GeV'^] 

xVdf 

^ (-^'max ) 

A(Z/ = oo) 

pow 

PS 

10,24 

0.570(35) 

41.2(7.1) 

31.65 

-5.57% 

-11.40% 

pow 

V 

10,24 

0.872(19) 

35.1(5.3) 

3.34 

-1.38% 

-4.87% 

pow 

N 

10,24 

1.255(43) 

88.5(9.3) 

5.83 

-2.99% 

-8.84% 

exp 

PS 

10,24 

0.624(13) 

65.9(4.2) 

5.59 

-2.41% 

-2.91% 

exp 

V 

10,24 

0.9125(92) 

63.5(5.6) 

1.19 

-0.17% 

-0.51% 

exp 

N 

10,24 

1.372(22) 

142.7(9.5) 

2.37 

0.18% 

-0.32% 


TABLE XII: Same as Table XI for (/I, k) = (5.6,0.158). 

Fit type 

H 

[Li,L2\ 

niH [GeV] 

c [GeV'*] 

xVdf 

^(Z/max) 

A(L = oo) 

pow 

PS 

12,24 

0.417(33) 

55.3(9.2) 

8.55 

-2.67% 

-14.81% 

pow 

V 

12,24 

0.780(59) 

62.5(13.1) 

5.22 

-2.10% 

-9.86% 

pow 

N 

12,24 

1.0894(92) 

124.0(2.3) 

0.07 

-0.44% 

-11.30% 

exp 

PS 

12,24 

0.466(20) 

47.9(4.2) 

3.86 

-1.50% 

-4.91% 

exp 

V 

12,24 

0.836(45) 

53.1(10.0) 

4.31 

-1.27% 

-3.41% 

exp 

N 

12,24 

1.208(20) 

104.1(5.1) 

0.50 

1.30% 

-1.65% 

TABLE XIII: Same as Table XI for (/3, k) = (5.32144, 0.1665). 

Fit type 

H 

[LuL2\ 

tuh [GeV] 

c [GeV'*] 

xVdf 

^(Z/max) 

A{L — (X)) 

pow 

PS 

12,16 

0.428(22) 

-15.0(16.3) 

2.09 

-0.83% 

2.23% 

pow 

V 

12,16 

0.743(32) 

23.3(29.0) 

0.77 

0.54% 

-2.08% 

pow 

N 

12,16 

1.037(63) 

90.9(52.8) 

1.87 

0.98% 

-6.05% 

exp {c = C) 

PS 

12,16 

0.4089(81) 


2.04 

-2.36% 

-2.36% 

exp 

V 

12,16 

0.756(20) 

18.4(26.7) 

0.86 

0.64% 

-0.31% 

exp 

N 

12,16 

1.088(42) 

74.9(50.5) 

2.32 

1.19% 

-1.46% 














